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ABSTRACT 



We report on an accurate measurement of the CXB in the 15-50 keV range 
performed with the Phoswich Detection System (PDS) instrument aboard the 
BeppoSAX satellite. We establish that the most likely CXB intensity level at 
its emission peak (26-28 keV) is ~40 keV cm -2 s _1 sr -1 , a value consistent 
with that derived from the best available CXB measurement obtained over 25 
years ago with the first High E nergy Astronomical Observatory satellite mission 
(HEAO-1; iGruber et al.lll999l ). whose intensity, lying well below the extrapo- 
lation of some lower energy measurements performed with focusing telescopes, 
was questioned in the recent years. We find that 90% of the acceptable solu- 
tions of our best fit model to the PDS data give a 20-50 keV CXB flux lower 
than 6.5 x 10~ 8 erg cm -2 s _1 sr _1 , which is 12% higher than that quoted by 
Gruber et al.l (119991 ) when we use our best calibration scale. This scale gives 
a 20-50 keV flux of the Crab Nebula of 9.22 x 10~ 9 erg cm~ 2 s _1 , which is in 
excellent agreement with the most recent Crab Ne bula measurements and 6% 
smaller than that assumed by I Gruber et al.l (Il999l ). In combination with the 
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CXB synthesis models we infer that about 25% of the intensity at ~ 30 keV 
arises from extremely obscured, Compton thick AGNs (absorbing column den- 
sity iV H > 10 24 cm -2 ), while a much larger population would be implied by 
the highest intensity estimates. We also infer a mass density of supermassive 
BHs of ~ 3 x 10 5 M Mpc~ 3 . The summed contribution of resolved sources 



( jMoretti et al.l 120031 ) in the 2-10 keV band exceeds our best fit CXB intensity 
extrapolated to lower energies, but it is within our upper limit, so that any signifi- 
cant contribution to the CXB from sources other than AGNs, such as star forming 
galaxies and diffuse Warm- Hot Intergalactic Medium (WHIM), is expected to be 
mainly confined below a few keV. 

Subject headings: X-rays: diffuse background — X-rays: general — galaxies: 
active — cosmology: diffuse radiation — cosmology: observations 



Introduction 



The cosmic X-ray background (CXB) is contributed mainly by active galactic nu- 
clei (AGN) powered by accreting supermassive black hole s at the centers of large galax- 
ies ( ISetti fc Woltjerl Il989l ; IComastri et al.l Il995l ; iGillil 12004 ) . Optically bright quasars and 
Seyfert galaxies dominate at low energies (up to a few keV), while obscured AGNs, w hich 
outnumber unobscured ones by a factor 3-4 fjUeda et al.l 120031 : lLa Franca et al.l 120051 ) , are 
responsible for the bulk of the CXB at high energies (>10 keV). However th e CXB intensity 



level is still a matter of debate. After the first pioneer CXB measurements ( jHorstman et al. 



19751 ). the major effort to get a reliable estimate of the spectrum in a broad energy band 
(2-400 keV) was performed in the late 1970's with the A2 and A4 instruments aboard the 
first High Energy Astronomical Observa tory (HEAO-1). The A2 results (3-45 keV) were 



first presented by iMarshall et al.l (119801 . hereafter M80), while the fi nal results ob t ained 



with the A4 Low Energy Detector (LED, 13-180 keV) were reported by I Grub er et al.l (11999 



hereafter G99), who also presented the conclusive results from both experiments. According 
to these authors the CXB energy spectrum J(E) in the 3-60 keV interval is well represented 
by a power-law (PL) with a high energy exponential cutoff (cutoffpl), while the corre- 
sponding E J(E) spectrum shows a characteristic bell shape with a maximum intensity of 
42.6 keV (cm 2 s srT 1 at 29.3 keV. 



After HEAO-1 there have been many other CXB measurements at low energies (<15 
keV) with both imaging and non imaging telescopes aboard satellite missions, but at high 
energies (>15 keV) no accurate measureme nt has been published yet. A major effort has 
been recently performed with INTEGRAL (jChurazov et al.l 120061 ). In this case the CXB 
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measurement was obtained by leaving the Earth disk (angular size of ~ 5.4° as seen by the 
satellite) to cross the much larger field of view (FOV) of the satellite mask telescopes and 
fitting the depression in the count rate with a multi-component model, in which one of the 
components was the CXB and the others were several, e.g., the Earth X-ray albedo, the 
instrumental background, the celestial sources in the FOV. The many assumptions about 
the contribution, at the epoch of the measurement, of all the components to the depression 
level and its dependence on energy and ti me, make this measurem ent somewhat difficult for 
an unbiased estimate of the CXB. Indeed IChurazov et al.l (120061 ) wish further INTEGRAL 
observations at other epochs "to verify the agreement of observations and predictions". 

In Table[TJwe compare the major CXB results. At low energies (1-15 keV), it is apparent 
a low spread of the PL photon index T and a high spread (up to ~40%) of the CXB inten- 
sity, exemplified by the ratio ^heao-T between the measured 2-10 keV intensities and that 
measured with HEAO-1, with the lowest CXB estimates obtaine d with HEAO-1 A 2 (M80) 
and the highest with the focusing telescopes aboard B eppoSAX (IVecchi et al.lll999l). XMM- 
Newt on ( Lumb et aPbood ; Ibe Luca fc Molendil 12004 ) and Chandra (jliickox fc Markevitch" 
20061 ). The r atio (= 1.15) obtained with the collimated Proportional Counter Array aboard 
Rossi-XTE flRevnivtsev et ajjjgOOgjt and that (= 1 .19) obtained from the re-analysis of the 
HEAO-1 A2 measurement (jRevnivtsev et al.ll2005l ) can be lowered to 1.04 and 1.07, respec- 
tivel y, by a more re liable calibration of the flux scale. Indeed the adopted 2-10 keV Crab 
flux (jZombeck!ll990l ) is ~11% higher than the mean value of all other Crab flux estimates. 
In the 20-50 keV energy band, which is common to the high energy X-ray experiments, the 
few available measurements show a low spread of the ratio -Rh°eao-i V but systematic errors 
in the CXB intensity estimates cannot be excluded. 



Driven by these discordant results several authors ( 
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la et al.ll2003: 


De Luca & Molendi 




2004; 


Comastri 


2004; 


Worslev et al.ll2005: 


Ballantvne et al. 


2006; 


Hopkins et al. 


2006; 


Worslev et al. 



20061 ). in their evaluation of the fraction of the CXB that can be resolved into individual 
sources or in their CXB source synthesis models, assume the 3-60 keV CXB spectrum ob- 
tained with HEAO-1 (G99) to be corrected in its shape, but underestimated in its intensity 
due to systematic errors in the absolute area calibration and/or instrumental background 
subtraction. As a result, the HEAO-1 CXB intensity is increased upward by a factor up to 
wl.3-1.4 over the entire energy band. 

In order to establish whether such CXB intensity renormalization is justified, and thus 
to constrain the size of a population of highly obscured AGNs and to infer the presence 



of ot her source populations and/or of a truly diffuse component (WHIM, see iKuntz et al. 



200ll ). we have performed an accurate measurement of the total (resolved plus unresolved) 
high energy (>15 keV) CXB intensity by exploiting the pointed observations performed with 
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the Phoswich Detection System (PDS) aboard the BeppoSAX satellite (IBoella et al.lll997al ). 



Instrument and calibrations 



A detailed description of the PD S instrument and its in-flight performance can be 
found elsewhere ( Frontera et al. 1997al lbh . Although the PDS was not designed to perform 



a measurement of the CXB (Field of View, FOV, of only 1?3 FWHM), a posteriori it was 
realized that such measurement would be possible given the very good performance of the 
instrument: high temporal stability, very low background level 5i 5 _ 3 oo kev = 1-6 x 10~ 4 counts 
(cm 2 s keV) -1 , high flux sensitivity and good energy calibration. The 15-300 keV limiting 
sensitivity corresponds to about 1% of the background level with a marginal influence (0.3%) 
of systematic errors in the background subtraction, also thanks to the continuous monitoring 
of the background with two rocking collimators which alternated with a default dwell time of 
96 s between the neutral position (ON-source) and two default symmetrical positions offset 
by ±3?5 (±OFF-source). 

The ON-axis PDS response function was determined by means of pre-launch calibra- 
tions combined with Monte Carlo calculations, and it was tested during the BeppoSAX life 
time (6 yrs) with 7 repeated observations of the Crab Nebula (twice in 1997 and 1999, and 
once in the years 1998, 2000, and 2001) that were performed simultaneously with the other 
Narrow Field Instruments (NFIs) on board. By assuming a power-law (pl) model of the form 
1(E) = NE~ r photons cm 2 s 1 keV 1 it was found that the photon index T obtained with 
the PDS remained unchanged in all of these observations within the statistical uncertainties, 
with a mean value of 2.121 ± 0.001 and a reduced x 2 /dof = 25.5 for 66 degrees of freedom 
(dof) in the 15-200 keV band. A better fit (x 2 /dof = 17.1 for 64 dof) in this energy band was 
obtained with a broken power-law (bknpl) with a mean low energy index T\ = 2.113±0.001, 
a break energy E\> = 74 ± 2 keV and a high energy index T2 = 2.198 ± 0.005. A bknpl 
spectrum was also the best fit to other broad band high energy measurements (see, e.g., 



Bartlettl Il994l ) . The still high x 2 /dof obtained with the bknpl model can be lowered and 
made compatible with the \ 2 statistics if a systematic error of only 1% in the used PDS 
response function is assumed. 



A cross-calibration of the PDS and MECS telescopes (see, e.g., IBoella et al.l Il997bl ) 
performed with the Crab has provided a time averaged normalization ratio at 1 keV between 
PDS and MECS of i?(PDS/MECS) = 0.928 ±0.001 and 0.917 ±0.001 under the assumption 
of a PL or a bknpl model, respectively. When corrected upward for this normalization ratio, 
the mean value of the PL normalization parameter is N = 9.54 ±0.01, while iV = 9.39 ±0.02 
for a bknpl. A comparison of this Crab spectrum with that obtained with other instruments 
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(see Fig. [T]) shows that this renormalized PDS spectrum is consistent within 8% with the 
extrapol ation of the 0.3-10 k eV spectrum obtained with the XMM-Newton EPIC-MOS 



came ra (jWillingale et al.l 1200 



19891 ) and with the classical 



) , with the Cra b spec trum obtained with HEAO-1 A4 (jJung 



Toor k, Sewardl (119741 ) results. In addition i t is in excellent 



agreement with the mean value of Crab spectrum quoted by iKirsch et al.l (120051 ). As can 
also be seen from Fi g. [fl with r espec t to the Crab 15-50 keV spectrum obtained with the 
PDS, that quoted by IZombeckl (Il990l ) is higher by ~ 30%, while that reported by iBartlett 
( 11994J ) is lower by ~ 15%. In the following we assume the renormalized Crab spectrum as 
our calibration scale for the CXB estimate. 

The OFF-axis response function of the PDS was tested with the Crab Nebula during 
the BeppoSAX Science Verification Phase. The Crab was observed in September 1996 at 
different offset angles with respect to the instrument axis and with a roll angle such as to get 
the narrowest a ngular response of the hexagonal collimators. Results of those measurements 
were reported (IFrontera et al.lll997al ) and now confirmed (see Fig. [2]). Up to 100 keV the 
angular response of each PDS unit to the Crab is well fit by a triangular function as expected. 
The derived FWHM is 9^ = 1?3 2, which is fully consistent with that derived in the pre- 
flight tests (IFrontera et al.lll997bl ). We also found that the response function is independent 
of the offset angle 9 apart from the exposed area through the collimators, which linearly 
decreases with 9. From the angular response of the PDS, we obtain (see Appendix A) for 
the geometric factor of the telescope a value G = 0.295 cm 2 sr, and for its solid angle a value 
Q = 4.624 x 10~ 4 sr. 



3. Measurement of the unresolved CXB 

3.1. Adopted method 

The measurement of the unresolved CXB photon spectrum stands as one of the most 
difficult tasks in observational X-ray astronomy. Among other things, it requires the knowl- 
edge, for each energy channel, of the instrument intrinsic background count rate v^ n to be 
subtracted from the total background v s ^ v measured during the observation of a blank sky 
field {y & ^ y = vcxb + ^in, where v C xb is the CXB count rate entering through the tele- 
scope FOV and v in is the intrinsic background). Systematic errors in the vqxb estimate can 
intervene either in a biased selection of blank fields or in a wrong z/j n evaluation or in both. 

In the case of the HEAO-1 A2 experiment, the total background level v s ^ v measured 
during the observation of a blank sky field was simultaneously observed through two colli- 
mators with different FOVs, one with a solid angle twice that of the other. Assuming U{ n 
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to be independent of the instrument FOV, the difference between the two total background 
levels removes the intrinsic bac kground and gives directly vcxb (M80). In the case of the 
HEAO-1 A4 experiment (G99. iKinzer et al.lll997l ). the unresolved CXB spectrum was de- 
rived by subtracting from the v s ^ y ', measured when the detectors observed blank sky fields, 
the background measured when the FOV of the detectors was shielded with a shutter made 
of a scintillator detector in anti-coincidence with the main detectors. In both cases some 
small bias in the vcxb estimate cannot be excluded given that the intrinsic background is 
dependent on the mass exposed to the environmental radiation: in the first case, there are 
two different collimator apertures, while, in the second case, the shutter could modify the 
intrinsic background. 

Our measurement of the unresolved vcxb count rate is based on the Sky-Earth Pointing 
(SEP) method, in which we subtract from the background level v s ^ v measured from a blank 



sky 



sky\ 

VCXB + V in ) 
.Earth 



the count rate level measured when the telescope is pointing 
v a + v fn Hh i where v& is the count rate due to the X-ray 



sky field 

to the dark Earth [y 

terrestrial albedo entering through the telescope FOV). The difference spectrum D(E) = 
{ycxB - va) + (v s in y - vf n aTth ) becomes D(E) = v C xb - va if vg arth = vf- Given that 
the radiation environment should not change looking to the dark Earth or to the sky, the 
latter condition is expected to be satisfied if both measurements are performed at the same 
cutoff magnetic rigidity and at the same time distance from the South Atlantic Geomagnetic 
Anomaly. The S EP strategy was als o adopted for the CXB measurement perfor med with 
the ASCA GIS fcushino et all I2OO2I ) and BeppoSAX MECS dVecchi et al.lll999h im aging 
telescopes, and with the Rossi-XTE PCA collimated detector (jRevnivtsev et al.l 120031 ) . 



3.2. Data selection 

In order to make sure that vf^ Tth = u^ v we performed a careful selection of the available 
data. The difference D(E) was derived only for Observation Periods (OPs) of > 10 ks 
duration during which the corresponding v s ^ v and i/^ arth were measured at similar mean 
values of the cutoff rigidity. In addition, given that the instrument mass distribution exposed 
to the sky (or Earth) and thus the intrinsic background can change with the collimator offset 
angle, we separately derived D(E) for the ON, +OFF and —OFF collimator positions. For 
the Earth pointings, we selected only those with the PDS axis well below the Earth limb. This 
method does not require a variable instrument configuration but requires the measurement 
of the albedo spectrum, which is not negligible at energies >15 keV. 

In order to satisfy the blank sky field condition, in addition to discarding all those 
pointings within 15° from the Galactic plane, for the OFF-source pointings we filtered out 
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those observations for which the +OFF and —OFF fields could be contaminated, e.g. from 
serendipitous X-ray sources, fast transients or solar flares. This selection was done by 
excluding from the sample those observations for which the difference between the two offset 
spectra (+OFF minus —OFF) were inconsistent with zero at 98% confidence level. For 
the ON-source pointings we accepted only those fields for which the difference between the 
ON-source count rate and count rate measured at either +OFF and —OFF is consistent 
with zero within la, and for which a fit with a null constant to the difference between the 
corresponding offset spectra (+OFF minus —OFF) gives a x 2 P er degree of freedom in the 
range 0.8-1.2. All the sky observations were done only when the instrument axis pointed at 
a direction at least 5° away from the Earth limb. 

As a result of the above selections, from the entire set of 868 BeppoSAX OPs off the 
Galactic plane, the number of useful OPs becomes 275 (127 ON-source, 71 +OFF-source, 
and 77 — OFF-source) with a total exposure time of 4031 ks. The dark Earth was observed 
for a total of 2056 ks. 



3.3. Results 

We obtained three difference spectra: D ON (E) for the ON-source pointings (2350 ks 
of exposure time), D +OFF (E) for the +OFF-source pointings (800 ks) and D_ OFF (E) for 
the —OFF-source pointings (881 ks). They are shown in Fig. [3j As can be seen, all are 
consistent with each other within their uncertainties, as also confirmed with a run test 



(see, e.g., iBendat fc Piersoll Il97ll ) . The small observed deviations from each other give an 
indication of the systematic errors made in the estimate of the difference v^ v — vfn Tth ■ 
However, the results from the offset and ON-source pointings agree within the statistical 
errors. On the basis of these results, for the derivation of the CXB intensity we used the 
sum D(E) = Don(E) + D + off(E) + D-o FF (E) which is well determined up to 50 keV (see 
Fig. ID . 

When D(E) is fit with a PL model, we find an unacceptable x 2 /dof (= 37.4/23) in the 
15-50 keV band, with a probability of 2.5% that this high x 2 value is due to chance. Instead, 
in the 20-50 keV band, an acceptable value (= 15.5/16) is found with a best fit PL photon 
index T at 90% confidence level of 3.2 ±0.2, much higher than that of the CXB (see Tabled]). 
Both results are expected if D(E) does give the difference spectrum between the CXB and 
the albedo from the dark Earth, with the first result mainly due to the presence of a low 
energy cutoff in the albedo spectrum (see Appendix B.l). The high T also unequivocally 
shows that the albedo spectrum, above its cutoff, is harder than that of the CXB. Thus, in 
the 15 to 50 keV band, we fit D(E) with the difference of two model spectra, one to describe 
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the unresolved CXB spectrum and the other to describe the albedo radiation spectrum. 

For the albedo model spectrum we used a photo-electrically absorbed power-law Ia{E) = 
exp (-t A fi A )N A (E/20)- rA photons cm 2 s 1 keV 1 sr 1 , where N A is the normalization 
constant at 20 keV, \i A is the air absorption coefficient in units of cm 2 /g, and t A is the 
atmospheric depth that describes the well known cutoff in the albedo spectrum at ~ 30 keV 
(see Appendix B.f ). To describe the albedo low-energy cutoff, we deve loped an atmospheric 



absorption model within the XSPEC software package (jArnaudlll996l ). that we adopted for 



our model fitting. The model makes use of a grid of values of air ma ss X-ray attenuation 



coefficients as a function of the photon energy (jHubbell Seltzerlll996l ). To model the CXB 



spectrum we assumed the CXB spectral shape obtained with HEAO-1 A2+A4 (G99). Thus 
we used as input models a cutoffpl and a PL which, in the 15-50 keV interval, still gives 
a good description of this shape. 

In order to better constrain the CXB spectral parameters, the fit of D(E) to the data 
was not performed by leaving free to vary all the models parameters. In return, different 
fits were performed, each one with different set of values of the parameters frozen in the 
fits. These parameters were: the CXB photon index T (when it was not left free to vary see 
Tab. [2]), the atmospheric depth t A and the albedo PL photon index T A . The allowed ranges of 
T were from 1.9 to 2.1 for the PL and from 1.2 to 1.4 interval for the cutoffpl, the allowed 
range of t A (in units of g cm -2 ) was from 1.4 to 6 and that of T A was from 1.3 to 2.0. Only 
the cutoff energy E c was always frozen to the value of 41.13 keV found with HEAO-1 (see 
Tabled]), given that in our energy band (15-50 keV) the fits were insensitive to E c . These 
ranges include the values obtained in past measurements and take into account the fact that 
the albedo mean slope is harder than that of the CXB, as previously discussed. In this way 
the entire parameter space was explored. By uniformly subdividing the allowed ranges of 
t A , T A and T in a certain number of subintervals, we obtained a grid of best fit values of 
those parameters that were left free to vary in the fits Q. From this grid we have derived 
the frequency distribution of one of the most important quantities that characterize our 
CXB estimate, i.e., the 20-50 keV integrated intensity Iq™#(20 — 50 keV) of the unresolved 
CXbB In Fig. [5] we show this distribution, for both a PL and cutoffpl CXB model. It was 
obtained by performing 1200 trials (20 steps for t A , 20 for T A and 3 for T), and accepting 
only those spectral solutions (532 in the case of the PL, 1091 in the case of the cutoffpl) for 
which the PL photon index of the reconstructed D(E) spectrum in the 20-50 keV band is in 



1 In all tables, those parameters that were frozen in the fits are shown in square brackets. 

2 We adopt this energy band given that it is common to all high energy X-ray experiments quoted in 
Table [T] and Tabled 
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the range 2.9-3.5, consistently with the observed slope of D(E). (Without this constraint the 
peak of the distribution occurs at slightly lower values of Iq^{20 — 50 keV).) We verified 
that the moments of this distribution are independent of the number of chosen steps. 

As can be seen, the frequency distribution of Iq^(20 — 50 keV) is well peaked (also the 
frequency distribution of the Ncxb and Na normalizations and of T show a similar shape), 
with a mean value slightly higher than that at the maximum of the distribution (see Fig. [5]). 
We have exploited this distribution to better constrain the range of the parame ter values that 



were frozen in the single fits. Having adopted the maximum likelihood method (jJanossylll965l ) 
for the best estimate of the model parameters, we considered as best fit parameter values, 
reported in Tabled those that, in the narrow flat top region of the frequency distribution 
around the mean value, give the minimum y 2 . In Fig. H]we show the best fitting curve to 
D(E) and the corresponding residuals in the case of a PL as input model for the unresolved 
CXB, and dark Earth albedo spectral parameters frozen at the values shown in Table [2j 

Using the J^xb s (20 — 50 keV) frequency distributions, we also derived the upper limit 
to the unresolved CXB estimate finding that, independently of the CXB input model, 90% 
of the data points have I^^{2Q — 50 keV) < 6.5 x 10~ 8 erg cm -2 s~ x sr _1 . We assume this 
value as upper limit to our estimate of the unresolved CXB. We notice that, by increasing the 
range of the CXB photon index values to 1.8-2.2 for the PL, and to 1.1-1.5 for cutoffpl, we 
find similar I^^{2{] — 50 keV) distributions with a mean at slightly lower intensity values 
but with upper limit almost unchanged with respect to that given above. 

In Appendix B.2 we report and discuss the obtained results on the terrestrial albedo 
from the dark Earth. For a comparison of the CXB spectrum with the albedo spectrum see 
Fig. El 



4. The total (resolved plus unresolved) CXB 

Exploiting the PDS pointings, we also performed an estimate of the contribution of 
resolved sources to the 15-50 keV integrated CXB intensity. This estimate could not be 
done using the ON-source pointings of the PDS, given that most of them were pointed 
observations of specific targets by the BeppoSAX Narrow Field Instruments LECS, MECS, 
HPGSPC, and PDS. However, 33 of the PDS OFF-source fields, that were excluded from the 
data set for the unresolved CXB intensity determination (see Sect. I3.2p . showed significant 
count excesses (2.5-7cr), consistent with the presence of serendipitous X-ray sources. The 
spectra of these excesses, when fit with a PL model, gave a weighted mean value of their 
photon indices equal to 1.65 ± 0.20, while their intensity gave a 15-100 keV energy flux in 



-10- 



the range (1.4-20) xlO -12 erg cm -2 s _1 . A search of their counterparts is now in progress. 
Their contributio n increases the CXB i ntensity by 4.7%. Adding the contribution of brighter 
sources (see, e.g. 



Krivonos et al.ll2005l ) does not significantly change this figure. 



Taking into account these results, in Table[3]we report, for each of the used input models, 
the normalization Nq^b °f the total CXB (unresolved plus resolved) spectrum, while in Fig. [6] 
we show, for the PL and cutoffpl models, the best fit E J(E) spectrum of the total CXB. 
In Table E] and Fig. O we also compare our measurement with the past results. 



5. Discussion 



It is apparent from Fig. [6] and Table [3] that our best fit Iq^ b (20 — 50 keV) is in excellent 
agreement with that obtained with HEAO-1 A2 (M80), and slightly lower (from 3 to 10%, 
depending on the input model) than that quoted by G99. (As discussed in Section [2J the 
use o f our fl ux scale calibration is more realistic; had we used the scale calibration given by 
June! Jl989h for HEAO-1 A4, our best fit /££ B (20 - 50 keV) would range from 0.95 to 1.03 
times the corresponding value derived from G99.) The best fit value of the maximum CXB 
flux density is obtained in the 26-28 keV band and ranges from 39.4 to 40.2 keV (cm 2 s sr) _1 
depending on the model assumed, with a statistical uncertainty in the centroid of ±1.5 
keV (cm 2 s sr)~ x at 90% confidence level for a single interesting parameter. 

We have also evaluated the upper limit to the CXB intensity that can be marginally 
accommodated by our data, by exploring the space of all the parameters involved in the fits. 
We find that, taking also into account the contribution of the resolved sources, independently 
of the CXB model, in 90% of this multi-parameter space the 7qxb(20 — 50 keV) is lower than 
6.8 x 10~ 8 erg cm -2 s _1 sr _1 , which is 12% higher than the best fit CXB intensity value quoted 
by G99 and 21% higher than that quoted by M80 (see Tabled]). 



Even this upper limit disagrees with the extrapolation to higher energies of the low 
energy (<10 keV) CXB estimates o btained with the focusing telescopes aboard BeppoSAX 



(IVecchi et al 



1999 ), XMM-Newt on (ILumb et al.ll2002l ; iDe Luca & Molendill2004h . and Chan- 



dra ( jHickox fc Markevitchl 120061 ). but it is consistent with those obtained with ASCA and 
RXTE. Thus, if we exclude a change in the CXB spectral shape derived with HEAO-1, our 
results raise the issue about the origin of the highest CXB intensities being quoted at lower 
energies. Differences in the flux scale calibration do not appear to be the origin of these 
discrepancies, as discussed in Section [2j One may think that part of the discrepant results 
could be due to the amount of sky solid angle surveyed, which is very large in the case 
of HEAO-1 and BeppoSAX PDS, and very small in the case of BeppoSAX MECS, XMM- 
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Newton and Chandra (see Fig. [7]), although, as also discussed by iBarcons et al.l (120001 ). this 
would imply that the sky regions surveyed by these telescopes are systematically brighter 
than the average sky sampled with HEAO-1 and BeppoSAX PDS. Another possible origin 
of the highest 2-10 keV CXB estimates could be due to systematic errors in the response 
function used for the diffuse emission (e.g., an underestimate of the stray light). For instance, 
in the case of MECS this function could be well tested and cross-calibrated with the PDS 
only for point-like sources. 

Independently of the CXB intensity issue at lower energies, our observational findings 
bear at least two important astrophysical consequences. Firstly, they provide a robust esti- 
mate of the accretion driven power integrated over cosmic time, including that produced by 
the most obscured AGNs. The AGN synthesis mo dels, tuned to attain the PDS CXB leve l 
and to account for the hard X-ray spectral shape (ILa Franca et al.ll2005l ; iGilli et al.l 120061 1 . 
predict that at the bright fluxes (> 10~ 12 erg cm -2 s _1 ) reachable by the PDS observa- 
tions and by coded mask instruments like INTEGRAL IBIS and Swift BAT, the detectable 
fraction of extremely obscured AGNs (Compton thick, Ah > 10 24 c m" 2 ) is ^10%, and i t 



increases to 20-25% at the fluxes reachable by focusing telescopes (e.g. lFerrando et al.1120 06'1. 



It sh o uld be noted that the highest CXB intensities c l aimed at low energies (IVecchi et al. 



Ferrando et al. 


2006). 


energies ( 


Vecchi et al. 



19991 : iDe Luca fc Molendil 12004 : iHickox & Markevitchl 12000 . see Table □]) w ould entail a 



much larger number (a factor 2-3) of Compton thick AGNs (IGilli et al.l 120061 ) . a prediction 
barely consistent with the present observational evidence. Our result implies a present black 
hole mass density of ~ 3 x 10 5 M Mpc -3 , using an admittedly uncertain bolo metric cor- 
rection of 30 for the 15-50 keV band and an efficiency of 0.1 in converting gravitational 
into radiation energy. This corresponds to a fraction of 6 x 10~ 5 of all baryons being locked 



into the supermassive black holes, using a cosmic baryon density of 
agreement with the "concordance" cosmology model. 



4 x 10 31 g cm 



m 



Secondly, under the assumption that the HEAO-1 spectral shape (G99) applies down 
to 2 keV, we fin d that the summed c ontribution of the observed X-ray source counts in the 
2-10 keV band (IMoretti et al.l 12003^ exceeds the PDS CXB best fit level by ~11%. This 
apparent contradiction vanishes if one takes into account the above discussed upper limit in 
the PDS CXB intens ity level and the error (±7%) associated with the source count evaluation 
(IMoretti et al.ll2003l ). As a consequence, our measurement suggests that it is quite possible 
that almost all the CXB in the 3-8 keV band has already been resolved into sources down 
to the faintest fluxes of the Chandra deep fields. Any substantial contribution to the CXB 
from other classes of sources and diffuse WHIM should be confined at photon e nergies below 



~3 k eV, as it has already been indicated in the case of star forming galaxies (IRanalli et al. 
20031 ) . 
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A. The solid angle of the PDS instrument 



The solid angle of the PDS is derived from the Geometric Factor of the telescope defined 
as (IPeterson et al.lll973l ; iHorstman et al.lll975l ) 



G = / A(8)dn, 



(Al) 



where A(6) is the exposed detector geometric area through the collimators at an offset angle 
9, and Q is the solid angle of the instrument Field of View (F OV). The expres s ion o f G, for 
an hexagonal collimator like that of the PDS, can be found in lHorstman et al.l (119751 ) and is 
given by 



G = nA(0) tan 2 (0.5tan 



_ x 0.939394 



h 



(A2) 



where A(0) is the ON-axis geometric area through the collimator (640 cm 2 ), d m is the 
diameter of the circumscribed circle to the hexagonal collimator cells, and h is collimator 
height. The ratio d m /h is related to the minimum FWHM of the angular response of the 
PDS 6*f w through the relation 

rl 9 

(A3) 



dm 2 



Using the value of 8f w derived from the offset Crab observations, we obtain a value of G 
0.295 cm 2 sr and, dividing by A(0) = 640 cm 2 , we find the telescope solid angle fl 
4.624 x 10~ 4 sr. 



B. The albedo spectrum 
B.l. The past measurements 



The terrestrial gamma-ray albedo radiation is mainly the result of the interactions with 
the upper atmosphere of the Cosmic Rays and, at X-ray energies, of the CXB and of the 
discrete X -ray source ra diation (Compton reflection). It was investigated since the late 



1960s (see |Petersonlll975l . for a review). Most of the properties of the atmospheric gamma- 



rays were obtained with balloon experiments, mainly launched from Palestine (Texas, USA). 
Empirical models of the atmospheric gamma-ray em i ssion, base d on observational results, 
were worked out by various authors (IPeterson! Il975l ; iLind Il975l ; iDean et al.l Il989l ). From 
these observations, it can be seen that the gamma-ray emission properties depend on var- 
ious parameters, like the geomagnetic latitude, the energy band, the direction of emission 
(downward, upward) and the altitude from the Earth. The atmospheric radiation measured 
by a satellite depends on these parameters, but, unlike the radiation observed by a balloon 
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experiment, a satellite observes only the radiation emerging from the top of the atmosphere 
(albedo radiation). Thus its spectral properties, specially at low energies (>10 keV), are 
expected to be different from those obs erved inside the atmosphere . Measurements of the 
X-ray albedo radiation are re ported in ISchwartz &: Peterson! (119741 ) for photon energies in 



the 10-300 keV band, and in Imhof et a. 



than 150 keV see, e.g., Letaw et al. (jl986 ) and Dean et al. ( 1989 ). 
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(j) above 4 keV , while for energies higher 



The albedo spectrum of ISchwartz fc Peterson! (119741 ). which was obtained with a 0.5 cm 
scintillator detector and a wide FOV (23° at zero response) aboard th e OS Q-3 satellite in a 
nearly circular orbit at an inclination of 33° and an altitude of 550 km (ISchwartz et al.lll9701 ). 
above 40 keV is consistent with a PL (see Fig. |SJ), while below 40 keV it shows a flattening 
with a definit ive low energy cutof f below 30 keV. This cutoff is also observed with balloon 
experiments (IPeterson et al.l 1 19731 ; iPetersonl Il975l ; ISchwartz fc Peterson! 119741 ) and it is at- 
tributed to self-absorption of the radiation collectively emitted from different atmospheric 
layers. The PL model above 40 keV is confirmed by the albedo spectrum measured with a 
50 cm 3 Ge(Li) cooled detector with a wide FOV as well (± 45° at zero respons e) aboard the 
low-altitude polar-orbiting satellite 1972-076B (see Fig. [8j llmhof et al.lll976l ). 



Following ISchwartz! (119691 ). the best fit to the >10 keV albedo spectrum measured with 
OSO-3is obtained with a photo-electrically absorbed PL Ia(E) = exp (— tAf^A)NA(E/20)~ rA pho- 
tons cm" 2 s _1 keV" 1 sr _1 , with photon index = 1.7 ± 0.3, and atmospheric thickness 
pA = 1.75 ± .15 g/ cm 2 (at 90% confidence level). In the case of the polar-orbiting satellite, 
Imhof et al.l (119761 ) found that above 40 keV the photon spectrum is consistent with a PL 
with index ranging from ~1.34 to ~1.39, depending on the latitude scanned. 

On the basis of these observations, we have assumed a photo-electrically absorbed PL 
as a model spectrum for the albedo radiation from the dark Earth. 



B.2. Our results 



In addition to the parameter values reported in Table El we show in Fig. [8] the derived 
spectrum of the terrestrial albedo from the dark Earth, compared with that of the CXB. 
As can be seen, the derived a lbedo spectrum is located between the OSO-3 results and 
the albedo spectrum derived by llmhof et al.l (119761 ). with a 20-50 keV integrated intensity of 



(8±2) x 10 -10 erg cm -2 s _1 sr _1 . It should be noticed that the intensity of the albedo radiation 
depends on the CR flux hitting the Earth, and thus on magnetic latitude. Given that in 
different Earth pointings we pointed to the Earth along generally different directions and 
thus to different magnetic latitudes, the derived spectrum is latitude averaged, as partially 
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done also in the case of the OS 0-3 and 1972-07 6B satellites due to their wide FOVs. We 
also notice that at different latitudes we do not observe the u pward albed o, but the albedo 
emerging at different zenith angles Z. However, as discussed by lLingl (119751 ). the atmospheric 
spectrum at low energies is expected to be not strongly anisotropic with Z. 
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Table 1. Summary of the CXB past measurements compared with HEAO-1 (G99). 



Instrument Ref. Energy band Model V kT or E c ^HEAo'-T ^HEAO-f 

(keV) (keV) 



Composite 


(1) 


1-20 


PL a 


1.59 ±0.02 




1.07 




Composite 


(1) 


20-200 


PL a 


2.040 ±0.013 






0.86 ±0.05 


Composite 


(2) 


20-165 


PL a 


2.17 ±0.07 






0.90 ±0.18 


HEAO-1/A2 


(3) 


3-50 


BREMSS 6 




[40] 


0.98 ±0.10 


0.92 ±0.05 


HEAO-1/A2 


(4) 


2-10 


PL a 


[1.4] 




1.19 ±0.06 




HEAO-1/A2 


(5) 


2-10 


PL a 


[1.558] c 




1.05 ±0.06 




HEAO-1/A2+A4 


(6,7) 


3-60 


CUTOFFPL d 


1.29 ±0.02 


41.13 ±0.62 


1 


1 


Rocket 


(8) 


2-6 


PL" 


[1.4] 




1.34 ±0.21 




ROSAT/PSPC 


(9) 


0.7-2.4 


PL a 


1.50 ±0.09 




1.20 ±0.05 




SAX/MECS 


(10) 


1-8 


PL a 


1.40 ± 0.04 




1.43 ±0.08 




ASCA/SIS 


(11) 


1-7 


PL" 


1.41 ±0.03 




1.06 ±0.05 




ASCA/GIS 


(12) 


1-10 


PL" 


[1.4] 




1.18 ±0.02 




XMM /EPIC-MOS/PN 


(13) 


2-8 


PL a 


1.42 ± 0.03 




1.30 ±0.14 




RXTE/PCA 


(14) 


3-20 


PL a 


1.42 ±0.02 




1.15 ±0.02 




XMM/EPIC-MOS 


(15) 


2-8 


PL a 


1.41 ±0.06 




1.36 ±0.10 




Chandra/ACIS-I 


(16) 


2-8 


PL a 


[1.4] 




1.33 ±0.13 





Note. — The reported energy band gives the interval in which the CXB spectrum and its model parameters have been 
determined. Rheao^l 

gives the ratio between the 2-10 keV intensity 1(2 — 10 keV) estimated from the reported parameters 

and that obtained with HEAO-1 (I H EAO-i(2 ~ lOkeV) = (5.41 ± 0.08) X 10~ 8 erg cm~ 2 s" 1 sr" 1 ). Likewise Rheao-i V S ives 
the ratio between the estimated 20-50 keV intensity 1(20 — 50 keV) and that obtained with HEAO-1(I HEAO _ 1 (20 - 50 keV) = 
(6.06±0.06) X 10~ 8 erg cm -2 s~ 1 sr" 1 ). In square parenthesis the parameter values that were kept constant by the quoted authors. 
Uncertainties are la errors for a single parameter. Where not reported, the uncertainties are very small or are not reported in the 
quoted papers. 



Refe renc es. — (11 [Horstman et alj 1 1975h; (2) iKinzer et ail dl978 j l; (3) M80; (41 IRevnivtsev et al.1 ll2005h; (5ljahoda fc et alj 
J2006h; f6)lGruberl ll992l l; ( 71 G99;78llMcCammon et al.l dl983ll; (91 iGeoreantopoulos et aljjl996h; (lO l lvecchi et al.l dl999h; (11) 



iGendreau et al.1 lll995ll; (12llKushino et al.1 j2002ll ; (13l lLumb et alj d200^ rr 141 IRevnivtsev et al.N20oi) ; (15llDe Luca fc Molend: ! 
J2004h ; (16l lHickox fe Markevitchl J2006h . 



"PL model: 1(E) tx £T r ph cm" 2 s" 1 keV" 



b Brcmsstrahlung as described in the XSPEC user manual where kT is the plasma temperat 
c The best fit was obtained with two PL plus a CUTOFFPL plus two edges. 
A I{E) oc E- r cxp (-E/E c ) photons cm" 2 s" 1 kcV" 1 sr" 1 



Table 2. Spectral parameters of the unresolved CXB and of the dark Earth albedo as derived from the PDS 

measurement. 



CXB model I Vggf c T E c d I A (20 keV) c T A t A e X 2 / dof /ggf (20 - 50 keV) 







Powcr- 


law a 




0.096 ± 0.003 


[1.98] 




0.014 ± 0.002 


[1.38] 


[1.91] 


9.43/23 


5 


62 ± 0.29 


Cutoflf 


power 


law 6 


0.151 ±0.005 


[1.4] 


[41.13] 


0.011 ±0.001 


[1.53] 


[2.42] 


9.2/23 


5 


27 ±0.29 


Cutoflf 


power 


law 6 


°- 160 -0.016 


1 4+U-4 


[41.13] 


011+ ,u,Ib 
u.uii_ 006 


[1.3] 


[1.4] 


9.0/22 


5 


62 ± 0.30 


Cutoflf 


power 


law 6 


0.141 ± 0.004 


[1.29] 


[41.13] 


0.0040 ± 0.0004 


[2.0] 


[4.98] 


9.2/23 




18 ±0.27 



Note. — ^cxb * s unres °l ve d CXB normalization at 20 keV, while I\ (20 keV) gives the intensity of the dark Earth albedo 
at 20 keV. The last column gives the unresolved 20-50 keV integrated intensity. The parameters that were frozen in the fits are 
shown in square brackets. The quoted uncertainties are errors at 90% confidence level for a single parameter. 

a / C XB(£) = V C XB(S/20)- r 

b /cx B (£) = N CXB (E/20)- F cxp(-E/E c ) 

to 

c In units of photons cm -2 s^ 1 keV^ 1 sr _1 ^ 

d In units of keV 

c In units of g cm -2 

'in units of 10 — 8 erg cm~ 2 s~ 1 sr~ 1 . 



Table 3. PDS results on the total (resolved plus unresolved) CXB compared to the past measurements. 



Experiment 


rtei. 


Energy band 
(keV) 


Model 


iV 


1 


kT or E c 
(keV) 


rtot /on c;n f 

7 cxb( 2U ~ 5U keV ) 


^Crab (2U — 50 KeV ) 


Composite 


(1) 


20-200 


PL a 


40.7 ±2.3 


2.040 ± 0.013 




5.2 ± 0.3 




Balloon 


(2) 


20-165 


PL a 


67 ± 13 


2.17 ±0.07 




5.47 ± 1.06 




HEAO-1/A2 


(3) 


3-50 


BREMSS 6 


13.95 ±0.70 




[40] 


5.6 ±0.3 




HEAO-1/A2+A4 


(4) 


3-60 


CUTOFFPL c 


7.877 ±0.08 


[1.29] 


[41.13] 


6.06 ± 0.06 


9.83 ±0.03 


SAX/PDS 


this paper 


15-50 


PL d 


0.100 ±0.002 


[1.98] 




5.89 ± 0.19 


9.22 ±0.01 


SAX/PDS 


this paper 


15-50 


CUTOFFPL e 


0.158 ±0.03 


[1.4] 


[41.13] 


5.52 ± 0.18 


9.22 ±0.01 


SAX/PDS 


this paper 


15-50 


CUTOFFPL 6 


0.167 ±0.017 


1.4 ±0.3 


[41.13] 


5.88 ± 0.19 


9.22 ±0.01 


SAX/PDS 


this paper 


15-50 


CUTOFFPL e 


0.148 ± 0.002 


[1.29] 


[41.13] 


5.43 ± 0.17 


9.22 ±0.01 



Note. — For each model, the energy band of the CXB spectral determination, the parameters of its photon spectrum and the 20-50 keV energy flux per steradian 
are reported. For flux scale calibration purposes, when available, also the 20—50 Crab flux predicted from the single experiments is reported. Uncertainties are la- 
errors. Parameters in square parenthesis are those kept fixed in the single fits. 

References. — (11 iHorstman et alj dl975h ; (2) lKinzer et al.l lll978h ; (3) M80; (4) G99 
*I(E) = NE~ r photons cm" 2 s^ 1 keV" 1 sr" 1 

b Brcmsstrahlung as described in the XSPEC user manual with N = K and kT is the plasma temperature. 
C I(E) = NE~ F exp (-E/E c ) photons cm" 2 s" 1 keV^ 1 sr" 1 
d I(E) = N(E/20)' r photons cm~ 2 s" 1 keV" 1 sr" 1 
C I(E) = N(E /20)~ r exp (-E / E c ) photons cm" 2 s" 1 keV -1 sr" 1 
In units of 10 — 8 erg cm -2 s~ 1 sr -1 
g In units of 10 — 9 erg cm -2 s _1 
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Fig. 1. — Upper panel: BeppoSAX 15-100 keV Crab spectrum (black points) compared 
with other measurem ents. Blue line: review of different measurement s in 5-70 keV 
jToor fc Sewardlll974h : Green line: GRIS balloon experiment in 20-1000 keV jBartlettlll994T ): 
Cyan line: IZombeckl (jl990l ); Purple line: collection of different measurements in 2-50 keV 
kirsch et al j|2005h; Oran ge line: 15-180 k eV HEAO-1 A4/LED measurement J.Tunglll989l ): 
Red line: XMM-Newton (jWillingale et al.ll200ll ). When not visible, error bars are smaller 
than the line thickness. Bottom panel: ratio between the Crab spectrum as measured by 
previous experiments and the PDS spectrum. With the exception of the Zombeck and the 
Bartlett measurements, all the other are consistent with each other within 8%. 
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15-100 keV PDS collimator response with Crab 
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Fig. 2. — Angular response of the PDS collimators for each phoswich unit measured from 
offset observations of the Crab in 15-100 keV. 
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Fig. 3. — Difference spectrum D obtained from the ON-source and OFF-source pointings. 
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Fig. 4. — Top panel: Average difference spectrum D = vcxb — v a of all the available data, 
along with one of the best fit models. In this case the CXB spectrum is modeled with a PL 
and the terrestrial albedo with an absorbed PL (see text). Bottom panel: residuals to the 
model. 
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Fig. 5. — Frequency distributions of the unresolved CXB 20-50 keV integrated intensity. 
The vertical dashed lines give the mean value of the distributions. Left panel: in the case of 
a PL CXB model with T in the range 1.9-2.1 (532 trials, see text). Right panel: in the case 
of a cutoffpl CXB model with T in the range 1.2-1.4 (1091 trials, see text). 
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Fig. 6. — Total (unresolved plus resolved) E J(E) CXB spectrum as observed with the PDS 
experiment (red points) compared with measurement results obtained with other missions. 
Upper panel: The energy spectrum J(E) modeled with a PL with T = 1.98 (see Table [2]). 
Bottom panel: The energy spectrum J(E) modeled with a CUTOFFPL with T = 1.4 (see 
Table EJ and E c fixed at the value of 41.13 keV obtained with HEAO-1 A2+A4 (G99). 
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Fig. 7. — Sky covered by the PDS compared with that covered by different focusing 
missions for the unresolved CXB determination. The 172 PDS poin tings correspond t o 
265 deg 2 . Sky c overage of focusing telescopes: 0.73 deg 2 (ME CS: IVecchi et al.l 119991 ): 



Kushino et al. 



2002); 



50 deg 2 (ASCA, 
deg 2 (Chandra; iHickox fc Markevitchl 120061 ) 



12 deg 2 (XMM-Newton; iLumb et al.l 120021 ): 0.5 
For the 34 XMM-Newton pointings used by 
De Luca fc Molendil (120041 ). coverage is 5.5 deg 2 , but coordinates are not available. 
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Fig. 8. — Comparison of the best fit photon spectrum of the dar k terrestrial albedo as derived 
by the PDS measurement (red points ) with the re s ults f ound by lSchwartz &: Peterson! (11974 
OSO-3 satellite, blue points) and by llmhof et al.l (119761 . 1972-076B satellite, green points). 
Also the best fit CXB photon spectrum assuming as input model a cutoffpl with T = 1.4 
(see Table E]) is shown for comparison with the derived albedo spectrum. 



